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We study elastic Na scattering and produce a quantitative correlation between the range of 
the effective potential and the energy of the system. This allows the identification of the waves and 
energies for which the scattering may be said to be peripheral. We then show that the corresponding 
phase shifts are sensitive to the tail of the NN potential, which is due to the exchange of two pions. 
However, the present uncertainties in the experimental phase shifts prevent the use of Na scattering 
to discriminate the existing models for the NN interaction. 
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I. INTRODUCTION 
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£NJ , In the present understanding of the NN interaction, long-range effects are ascribed to single pion exchanges, 
intermediate components are associated with exchanges of two and three pions, p's and u>, as well as baryon excitations, 
such as the A. At short distances, quark dynamics is supposed to dominate. 

The long-range one-pion exchange potential (OPEP) became well established in the 1960's as a necessary 
component of all realistic forces. This achievement was followed by an effort to determine the next layer of the 
interaction, associated with the two-pion exchange potential (TPEP), and various strategies were proposed to describe 
it. Early in 1971, Brown and Durso [0 pointed out that this component of the force is directly related to off-shell 
ttN scattering and hence strongly influenced by chiral symmetry. In a subsequent paper, Chemtob, Durso, and Riska 
H explored this relationship and investigated the phenomenological features of the TPEP obtained by means of 
dispersion relations. A very important step in this research program was the construction of the Paris potential [0, 
. where dispersion relations were used to relate empirical information about the irN process to the intermediate part 
of the force. This potential is rather successful in describing experimental data. 

The intermediate part of the potential may also be studied in the framework of field theory. In this case, one 
writes down a Lagrangian involving the relevant degrees of freedom and then evaluate a certain number of Feynman 
diagrams. This leads to an amplitude which is afterwards transformed into a potential. An important early work along 
this line was that of Partovi and Lomon ||, who considered a Lagrangian containing just pions and nucleons with 
a pseudoscalar (PS) coupling and evaluated the box and crossed box diagrams contributing to the exchange of two 
uncorrelated pions. A detailed study of the same diagrams using a pseudovector (PV) coupling was performed later 
by Zuilhof and Tjon H. The inclusion of the exchanges of resonances and baryon excitations led to the construction 
of the Bonn potcntial^D > which is also quite successful in reproducing experimental data. 

In a parallel line of development, several phenomenological potentials were constructed, which are also able to 
account for experimental information with the aid of adjusted parameters ]9]-|l2|]. 

An interesting feature of all approaches to the intermediate part of the potential is that their theoretical reliability 
decreases as one moves from the outer to the inner region. The intermediate range interactions exhibit a marked 
spatial hierarchy. The spatial features of a given interaction are determined by the mass exchanged in the t channel. 
In the intermediate part of the potential, the lightest system that can be exchanged involves just two pions and has a 
mass of about 300 MeV. Other important effects, associated with resonances such as the p and the lj are short ranged, 
because these states have masses around 750 MeV. When field theory is used, predictions for the inner parts of the 
mesonic sector of the potential are heavily influenced by form factors, bringing a lot of uncertainty to calculations. 
In the case of dispersion relations, on the other hand, predictions for the inner region are based on data which need 
to be extrapolated far away from their experimental region. 

Recently there has been a renewal in the interest on the intermediate part of the NN interaction, motivated by 
the realization that chiral symmetry provides a suitable theoretical framework for the calculation of a strong process 
|p3|-|l6|. In its minimal version, the chiral intermediate potential is based on a system containing just pions and 
nucleons. The theoretical foundations of this part of the potential are rather well established and it is reasonable 
to expect that it should become a standard ingredient of any modern NN potential. A shortcoming of the minimal 
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chiral potential is that it fails to reproduce experimental information in the case of the intermediate irN amplitude 
entering in the TPEP. In order to overcome this difficulty, one may extend the chiral model, so as to encompass other 
degrees of freedom. This possibility was recently considered by Ordonez, Ray, and van Kolck E^], who have shown 
that the inclusion of A's in the model improves its predictive power. Alternatively, one may choose to introduce the 
empirical information that is missing in the intermediate TtN amplitude in a model independent way, with the help 
of the Hohler- Jacob-Strauss [[L8| subthreshold coefficients [ ^9|j20| l. This led to a TPEP which yields a satisfactory 
description of NN data for waves and energies associated with peripheral scattering 21 22|. 

In spite of the considerable amount of activity related to the intermediate part of the NN potential in the last twenty- 
five years, no consensus was reached about the fine details of this component of the force. The dynamical content 
of the various models is not uniform and the profile functions they yield for different components of the potential 
do not agree. This picture poses the problem of defining criteria for establishing the merits and shortcomings of the 
various existing potentials, so that a choice can be produced. An obvious criterion of choice is the ability a given 
potential has of reproducing experimental data. This possibility has been on for a long time and does not work as a 
distinguishing criterion, because all existing modern potentials are able to explain well experimental NN phase shifts, 
but do so with the help of several free parameters that are adjusted ad hoc. Another problem about relying on NN 
observables is that the OPEP contributes to all channels and waves, making it difficult to isolate unambiguously the 
contribution of intermediate range dynamics. It would therefore be interesting to find ways of obtaining information 
about the intermediate part of the potential directly from empirical data. 

In this work we speculate about the possibility of obtaining such an information from the study of nuclcon-alpha 
(Na) scattering. The a particle is a rather suitable system for the study of the intermediate part of the nuclear 
interaction, because it is a boson that has no spin and isospin and hence cannot couple to a single pion. This means 
that the outer part of the nuclear potential surrounding the a involves two uncorrelated pions. The scalar-isoscalar 
channel of the two-pion exchange NN potential has strong central and spin-orbit components, and the same happens 
with the Na effective interaction. This picture has empirical support, since JVq experimental phase shifts show 
that central and spin-orbit effects are very important [^3 26 2^j24|]. Therefore we may expect that the peripheral 
scattering of nucleons by the a should be heavily dominated by two-pion exchange process and hence reflect the 
various approaches adopted in different NN models. When low-energy protons are used as probes, the Coulomb 
barrier prevents short distance interactions. Thus the effects due to the nuclear force show themselves as deviations 
from Coulomb amplitudes, similarly to what happens in low-energy nucleus-nucleus collisions p7| , p8[ . For instance, 
in sub-Coulomb Pb+Pb Mott scattering the intermediate range NN interaction accounts for deviations observed 
experimentally p7| ]. In the case of low-energy neutrons, on the other hand, angular momentum may be used to select 
the various regions of the potential. 

This work is organized as follows. The effective JVa potential is defined in sect. O and explicitly constructed in 



sect. III. In sect. IV we display our dynamical equations, which are based on the variable phase method. Finally, in 
sect. V] we discuss our results. 



II. DYNAMICS 



In this section we present briefly the main equations used in this work. In a complete treatment of the Na scattering 
problem, one has to deal with a rather complex five-body system. However, here we are interested in estimating the 
effects of the tail of the JVJV interaction, which manifests themselves at large distances and high values of the angular 
momentum. We therefore assume that the a subsystem remains undisturbed during the interaction. 

The wave function for the four-body nuclear system 4 He can be written, in terms of the spatial, spin, and isospin 
degrees of freedom, as 

\N% ■ ■ ■ Ni) — \n ■ ■ ■ r 4 ) (g) |spin) ® |isospin) 

= |R)|a), (2.1) 

where the collective variable R is given in terms of the individual coordinates r^: 

i 4 

i=l 

and | a) represents the 4 He ground state. This wave function is known to have T — 0, J 71 " = + , and both S and D 
spatial components [^9|. In order to simplify the calculation we neglect D waves, adopt a Gaussian structure for the 
spatial wave function |2S|] and write 
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= 7V Q exp ^ 
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(2.3) 



where = ( — r^), N a is a normalization constant, a is a parameter extracted from Ref. 29 1, and \xa) is the 
spin-isospin wave function. 

The dynamics of the four body system is determined by the Hamiltonian 



^ 2m ^ 

i—l 



V 2 

2M a 



(2.4) 



where 0^ = kr?/2 is the harmonic potential and k — 4a 2 /m. 

In order to isolate the motion of the center of mass, it is useful to use the Jacobi variables, given generically by 



Pi =53 a « r J' 



where the coefficients ay are summarized in the matrix 



a = 



1/2 -072 
1/6 0/6 -y/2/3 



0/12 0712 0/12 -0/4 



In terms of these new coordinates, the Hamiltonian for the isolated four-body system is 

Hi 



2M a 

where — Vf l /2M a describes the center-of-mass motion and 



V 2 



2m,- 



2kp\ 



is the intrinsic four-body Hamiltonian. The Schrodinger equation for the a particle is 

H a | a) = E a | a) , 

where E a = M a — 4m and |a) is now given by 



/ 4a \ 9/4 { * 2 V . 
| a) = cxp \-2a^p l J \x a ). 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



(2-10) 



The strong interaction between the incoming particle, hereafter labeled by o, and the nucleon i within the a particle 
is denoted by V i(r oi ). The study of the effects of this interaction is the main object of this work. 
The Schrodinger equation for this five nucleon system is 



H 5 \N ---N 4 )=E\N ---N i ), 



(2-11) 



where 



i=0 



= H 



V 2 
2m 



i>j — 1 i—l 



2M a 2m 



(2-12) 
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Approximating the five-body wave function by a product of the different clusters constituting the system 

\N ---N 4 )~\N ,TL)\a), 

we can rewrite the dynamical equations as 



H 5 \N ■ ■ ■ N 4 ) 



V 2 V 2 4 



\N ,R)\a) 



Introducing two new Jacobi variables 



= E\N ,TL)\a). 



s= -(r + 4R), 
x = r G - R, 



we have 



2(M a +m) 2[i j-^ 



V 2 4 



\N ,TL}\a) = (f?-£7 a )|JV ,R)|a), 



where 



M a m 4 
M = tt — : — —m. 



M a + m 5 

Writing down explicitly the projectile wave function as 

\N ) = \r )\ Xo ) 

and going to the CM. frame of the five-body system, we obtain 



V 2 4 



|x)| Xo )|a)=i; x |x)|x )|a), 



where E x = E — E a . Multiplying this equation by (a\ and integrating over the a coordinates, we get 

|x)|Xo) = E x \x)\ Xo ), 



V 2 



where 



W(x) = (a\^v oi ( *oi) I a) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



is the effective N— a potential. 



III. THE EFFECTIVE POTENTIAL 



The two-body strong potential V Q i is written as 



Voi = J2 V o T i S ( r oi)Pj i P t 



S 
oi ' 



S,T 



(3-1) 



where the indices T and S represent the total isospin and spin of the NN system whereas and are the 
corresponding projection operators. The explicit forms of the isospin operators are 



4 



P 1 



h_ T (o). T Wl 



3 + T (o) . T (i) 



(3.2) 
(3.3) 



and similar expressions hold for the spin degrees of freedom. 

The various radial components V^ s have the following structure 



V™ = VP{v oi ) + V[i(r oi )L 

+VP(v ol ) (3aW ■ r oi <rW • r oi - ■ «r«) 

+small components. (3-4) 

In this expression, L ^ is the orbital angular momentum for the pair oi, whereas V^ s , V[g , and Vrf s represent the 
central, spin-orbit, and tensor components of the potential, respectively. 

In the evaluation of the effective potential, one needs the expectation value of V i between the wave functions of 
the a particle. Noting that (x a |f^|Xa) = (Xa| T ^IXa) = 0; we have 



(Xa\V i(r i)\x a ) = Vc(r i) + V L s(r oi )L 



T {o) 

IT' 



where 



and 



V c (r oi ) = - [VS°(r oi ) + 3^°(r oi ) + ZV°\r oi ) + QV^ 1 ^)] 



V L s{v ol ) = -[Vll{v m ) + WlUr m )}. 



(3.5) 
(3.6) 
(3.7) 



The strong effective potential is then given by 

4 / , n 9/2 



W{x) = y—J J d Pi d P2 d P 3 ex P ( -2" Po 



Vc(roi) + V LS (r oi )L c 



A<>) 



cxp -2aJ2pj 



3=1 



In the evaluation of this expression, it is convenient to note that 

3 

r t = R + 4 



T 
2 ijPj 



3 = 1 



and hence 



v oi =x-J2aJjPj 

3 = 1 



(3i 



(3.9) 



(3.10) 



The angular momentum operator L oi is given by L oi = r oi x p oi , where p oi = —iV i- The use of the chain rule 
allows us to write 



L„, = x 



Y a fjpjj x ^Px-E a ^j 
i i 3 3 

J XXPx_ lE a ij(Pj X Px ) ~ X P*) 



3 = 1 

3 3 

mi a % a ^k(p3 X P*)- 

fe=l j=l 



k=l 



(3.H) 
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When the operator acts on the spatial wave function, we have 

Pfcexp ^-2«X]p^ = -*(- 4 «)Pfc ex P (^ 2a Yl p i 
Defining L x = x x p x , we write 

f <r(°) 
IT(x) = W c {x) + -W LS (x)L x • — . 



(3.12) 



(3.13) 



The Gaussian structure of the wave function allows several integrations to be performed analytically and we obtain 

1/2 



(16ct \ 1 f \ 
—J -J duuV c {u) |exp 



16a. 



exp 



16a, 
— (x + u) 2 



and 



/16a\^^ 1 f°° 
W LS (x) = 4 ( — \ -2 J duuV LS {u) x 



(3.14) 



(3.15) 



32ax 



exp 



16a. , 9 
— {x~u) 2 



32ax 



exp 



16a, , 
—(x + u) 2 



For a channel with total angular momentum J and orbital angular momentum L, we obtain the following form for 
the effective potential: 



Wj L (x) = W c (x) + ^ 



J(J+l)-L(L + \)-- 



W LS (x). 



Thus we have 



w l-Il( x ) = Wc(x) - -(L + l)W LS (x), 



W L+hL (x) = W c (x) + -LW LS (x). 



(3.16) 

(3.17) 
(3.18) 



IV. THE VARIABLE PHASE APPROACH 



We are interested in the properties of the tail of the NN interaction and hence it is convenient to calculate the 
Na phase shifts by means of the so-called variable phase method, as described by Calogero |5(| . This method is fully 
equivalent to the Schrodinger formalism, and has the advantage of providing a clear picture of the spatial influence 
of the potential over the phase shift and of its relationship with the centrifugal barrier. For the reader's convenience, 
we hereafter collect the main formulas needed in this calculation. As we disregard tensor interactions, there are 
only uncoupled channels. The wave ujl{t), representing a system with total and orbital angular momenta J and L, 
respectively, is written in terms of the radial Green's functions as 

ujl{t) = c JL {k,r)j L (kr) - s JL (k,r)h L (kr), (4.1) 

where k = \/2[iE x , ji(fcr) and fii(fcr) are the usual spherical Bessel and Neumann functions multiplied by kr whereas 
cjL{k,r) and sji,(k,r) are defined in terms of Vjl{t) as follows: 

c JL (k,r) = l-\ dr'Vj L (r')j L (kr')uj L (r'), (4.2) 

sj L (k,r)=-^-J^ dr'Vj L {r')n L (kr')u,, L {r'). (4.3) 
The phase shift Sjl (k) is related to these functions by 
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tan 5jl (k) = lim 



sjL(k,r) 



cjh{k,r) 

This motivates the definition of a variable phase Sjl (k,r) through 

sjL(k,r) 



Sjl (k, r) = arctan 



cjh{k,r) 



(4.4) 



(4.5) 



This function vanishes at the origin |p0|] and becomes the phase shift when r approaches in finit y. It is determined 
dyna mically by means of a first-order differential equation, obtained by differentiating Eq. (4.5), using Eqs. ( |4.2| ), 
(4.3), and manipulating the result to obtain 



^5 JL (k,r) = -\w JL {r)P 2 L [kr,S JL (k,r)} , 

where Pl is the uncoupled structure function, given by 

Pl [kr, Sj L (k, r)] = cos [S JL (k, r)] Jz(fer) - sin [S JL (k, r)] n L (kr). 
The integral expression for the variable phase is 



Sjl (k,r) 



1 

k Jo 



dr'Vj L [r')Pl\kr',8 JL {k,r')\ 



For future use, we define the ratio pji,(k, r) as 



PJL{k,r) 



Sjl (k,r) 
Sjl (fc,oo) ' 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



It represents the fraction of the phase shift that is generated by the part of the potential between the origin and the 
point r. 



V. RESULTS AND DISCUSSION 



The purpose of this work is to discuss the possibility of using Na scattering to study the tail of the two-pion 
exchange nucleon-nucleon potential. As stressed previously, the a particle is a spin-isospin scalar and hence it is very 
suited to the study of the intermediate-range part of the nuclear interaction. We therefore consider several realistic 
NN potentials g@,gjl0|^|@ and fold them into the a wave function, in order to obtain effective interactions, 
which are then used in the Schrodinger equation. The effective Na potential is dominated by scalar and isoscalar 
exchanges, giving rise to central and spin-orbit contributions, displayed in Fig. 13. In all cases, as expected, the 
spin-orbit component falls faste r than the central one at large distances. Both components are negative there and 
hence, according to Eqs. (3.17) and (3.18), they interfere destructively when J = L — i and constructively when 
J = L + ~. This means that phase shifts for peripheral waves of L are positive and larger for the latter class of waves. 

We arc interested in effects of the tail of the potential and hence it is important to establish, for each wave, a 
correlation between the energy and the region of the potential that determines the phase shifts. With that goal in 
mind, we have employed the variable phase method to evaluate the ratio pjL{k,r), defined in Eq. (4.9), for various 
energies, in order to determine the radii Re, and i?io for which pj^^k, r) is always less than 5 and 10 %, respectively. 
This means that 95% of the phase shift is generated in the region where r > R5 and, similarly, 90% for r > Riq. 
Results for the Argonne potential |l2| are presented in Table |[ where it is possible to notice some interesting features. 
The first of them is that, for a given value of L, the phases are slightly more sensitive to the short-range part of the 
potential for the channel with J = L — i. This feature of the problem can be fully appreciated in Fig. |2|, that shows 
the ratios pf 5/2 an d Pf 7/2 a t 15 MeV, which are representative of the other peripheral waves. The behavior of these 
curves are determined by the centrifugal barrier, responsible for the plateau close to the origin, combined with an 
effective potential which is repulsive at short distances and then attractive. 

The a particle has a rms of about 1.6 fm and, in a conservative approach, we may define the tail of the effective 
potential as beginning at 2.0 fm. Table || then informs us that such a tail of the Na potential determines more than 
95% of the phase shifts 8q 7 Sc g/2 up to 50 MeV and Sh 9/2 , Sh 11/2 up to more than 80 MeV 



In Fig. H we display the phase shifts for waves with 2 < L < 5. Inspecting the figures corresponding to the waves 
F, G, and H, one learns that the various potentials provide a coherent pattern, since all curves have the same general 
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trends. The phases are sensitive to the tail of the effective potential and hence one could, in principle, discriminate 
the medium range features of the NN interactions. In Fig. |3| we also include the available experimental values, taken 
|. The picture provided by these data, claimed to be precise, does not allow conclusions, for 
there are big discrepancies among phase shift analyses made by the various authors. So a drastic improvement in the 
precision of Na phase shifts is needed before they could be used to discriminate NN potentials. 

As a final comment, we would like to stress that our calculation is a very simple one and a more detailed study 
is required on the theoretical side. One of the points that deserves further attention concerns the construction of 
the a wave function, that should include D components and, if possible, be calculated using realistic NN potentials. 
Although it is reasonable to assume that intermediate states involving the 3 Li5, 2 Hes, or a excitations contribute 
little to peripheral scattering, it is also important to clarify this aspect of the problem. Finally, it would also be 
important to establish the quantitative role of the two-pion exchange three-body force, since its range is comparable 
to that of the TPEP. We are now investigating some of these questions and prefer to take the results of the present 
work as indicating reliably only that Na scattering allows the discrimination of the medium-range content of NN 
potentials. In order to be able to select one or a class of potentials as being better than others one needs both more 
precise calculations and experimental phase shifts. 
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3.26 


3.54 


20 


1.74 


2.02 


1.86 


3.02 


3.10 


3.26 


2.94 


3.12 


3.30 


3.56 


3.18 


3.44 


30 


1.56 


1.78 


1.68 


1.98 


3.06 


3.20 


2.88 


3.02 


3.20 


3.40 


3.06 


3.28 


40 


1.40 


1.58 


1.52 


1.74 


2.14 


3.14 


2.82 


2.96 


3.12 


3.30 


2.96 


3.16 


50 


1.26 


1.40 


2.23 


2.95 


1.94 


2.34 


2.80 


2.92 


3.06 


3.22 


2.90 


3.08 


60 


1.08 


1.20 


1.97 


2.91 


1.80 


2.06 


2.00 


2.88 


3.02 


3.16 


2.84 


3.00 


70 


0.78 


0.86 


1.12 


1.26 


1.68 


1.88 


1.84 


2.86 


3.00 


3.12 


2.80 


2.94 


80 


0.86 


1.06 


0.98 


1.08 


1.56 


1.74 


1.72 


1.98 


2.26 


3.08 


2.76 


2.90 



TABLE I. Radii -R5 and _Rio, in fm, as function of the energy E, in MeV, for which the normalized phase shift, calculated 
for the Argonne potential, is always less than 5% and 10%, respectively. It is worth noting the non monotonic behavior of the 
F5/2 wave, which is due to a change in sign in the Argonne phase shift just before 50 MeV. 
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FIG. 1. Effective Na central (continuous line) and spin-orbit (dashed line) potentials for the Paris [Q, Bonn (OBEPR) Q, 
supersoft core C Q, dTRS |10| , Argonne [ ]12| , and chiral interactions. 



FIG. 2. Phase shifts normalized to 1 [Eq. 4.9 at 15 MeV, as function of r for the waves -F5/2 (dashed line) and Ft/ 2 
(continuous line). 

FIG. 3. Na phase shifts for the Paris (P), Bonn (B), supersoft core C (S), dTRS (TV L Argonne (A), and chiral (C) AW 
interactions. The experimental results are taken from Refs. j§^](o) , j24| (o) , J2f|(+), and |2q1 ( x ) . 
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